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We show that universal functions play an important role in the observation of the 
fractal structure of space-time in the numerical simulation of quantum gravity. 



1 Introduction 

The fractal structure of space-time is one of the most interesting aspects in 
quantum gravity. Since we know well-defined theories of quantum gravity only 
in two dimensions, we here confine ourselves to the two-dimensional quantum 
gravity. We study the fractal structure of space-time from the points of view 
of the geodesic distance as well as the "time" in the diffusion equation. 
The partition function for the two-dimensional quantum gravity is 

V[9am 5{[v9-V)e- S ^ I**, (1) 



Vol(Diff) 

where <fi symbolizes matter fields. We have here fixed the volume of space-time 
V = fyfg by 8- function. The usual cosmological term A-fy/g is recovered by 
the Laplace transformation: Z\ = J Q dVe~~ AV Zy . 



1.1 Geodesic distance 

We first study the fractal structure by using the geodesic distance. Let us 
consider the total length of the boundaries Sv(R), which are separated from 

"This article is based on rcfs. jjjj and [p| which have been done in collaboration with J. 
Ambj0rn and J. Jurkiewicz. 
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a given point with a geodesic distance R. The average of Sv(R) is given by 

(Sv(R)) = (2) 

where Gy{R) is the two-point function defined by |l|] 



dg(£,£,o) denotes the geodesic distance between two points labeled by £ and £o- 
Since dR (Sv(R)} is a unit volume, we have 

JO 

dR(S v (R)) — V . (4) 
According to the scaling arguments, the following generic behavior is derived: 
(Sv(R)) = V X - U U{X) with X = ^. (5) 

Here a parameter v is introduced so as to let X = R/V u be dimensionless. 
Expanding (^) around R ~ 0, a kind of fractal structure is expected, i.e., 

(S V (R)) - const. v l - vdh R dh - x for R - , (6) 

which defines the intrinsic Hausdorff dimension dh ■ For a smooth d-dimensional 
manifold we have 1/v = dh = d. If (Sv=oo(R)) is non-zero finite, i.e., 

vd h = l, (7) 

a kind of "flat" fractal space is realized when V — > oo|] We call the fractal 
which becomes a "flat" fractal space for V — > oo by the "smooth" fractal| 

We will show later that the function J7(X) introduced in (|5|) plays an 
important role in numerical simulation. From (|5|) and (^|) we obtain 

U{X) = V v - l {S v {R)) = X^FiX), (8) 



'A counter-example exists. We find that m-th multicritical branched polymer model ir 
does not satisfy the condition (Q). 

c The fractal structure of space-time with infinite volume is discussed in |4 1. 



where F(0) is non-zero finite. From we find that U(X) is normalized as 

dXU(X) = 1. (9) 



In the case of pure quantum gravity, one can calculate Gy{R) analytically 
|jj using the transfer matrix formalism Q: (See also the appendix in ref. 0) 

H R) = l^Ar^e^^g^, ,10) 

while Zy = A/V 7 ^ 2 , where A is constant. Using (|J) and (|l0|), we find 

(SV(-R)) - const. R 3 for J? ~ . (11) 

Thus, we find X/u = dh = 4 in pure gravity. Note that the "smooth" fractal 
condition (^|) is satisfied in this case. 

1.2 Time of Diffusion 

Now, let us introduce the "time" of diffusion in order to analyze another aspect 
of fractal structure. The diffusion equation has the form, 



d_ 

dT 



K g %faT) = A g K g (Z,Z ;T). (12) 



We here consider the following initial condition, 

K g (Z,t 0] 0) = J(£-e„). (13) 

K g (£, £ ; T) is the probability of diffused matter per unit volume at coordinate 
£ at time T. We introduce the following wave function Ky(R; T), which is the 
average of K g (£,£ ;T) with a fixed geodesic distance R: 



Since Ky(R]T) satisfies 

/>GO 

/ dR(S v (R))K v (R;T) = 1, (15) 
Jo 
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the scaling arguments lead to 

K V (R-T) = yP(X-Y) with X = A y = 1 (16) 

where a parameter A is introduced so as to let Y = T/V x be dimensionless. 
Similarly to <M), we find 



/>oo 

/ dXU{X)P{X-Y) = 1 
Jo 



(17) 



Thus, P(X; Y) as well as U(X) are considered as a kind of universal functions. 

Now we consider Kyifi', T) which is the average of the return probability 
of diffused matter at time T. Around T ~ 0, K v (0; T) will be 

y\d B /2-l 

K v (0;T) ~ const. ^ for T~0, (18) 

which defines the spectral dimension d s . We also consider the average geodesic 
distance travel by diffusion over time T, 

(R V (T)) d i f f dR{S v {R))RK v {R-T) ~ const. V u - Xa T° for T~0. 
Jo 

(19) 

For a smooth d-dimensional manifold we have 2/X = d s = d and a = 1/2. If 
the surface is a "smooth" fractal, the right-hand sides of ([l8]) and ([ii]) take 
non-zero finite values for V — ► oo. Then we find 

Ad s = 2 , Act = 1/ , (20) 

as "smooth" fractal conditions. The universal functions for (|l8| ) and ( |l9|) are 
P(0;Y) = F^y(0;T) and J °° dXt/(X)XP(X; F) = V"-" (Rv(T)), respec- 
tively. 



2 Numerical Simulation 

The setup of the numerical simulations in ref M is as follows: We use the dy- 
namical triangulation, i.e., ensembles of surfaces built of equilateral triangles 
with spherical topology and a fixed number of triangles. These surfaces can 
be viewed as dual to planar (f> 3 diagrams. In the language of the <f> 3 theory, 
we include diagrams containing tadpole and self-energy subdiagrams. In sim- 
ulation both the standard flips and the new global moves (minbu surgery) are 
organized. The new move helps to reduce the correlation times We 
study various sizes of systems with 1000, 2000, 4000, 8000, 16000 and 32000 
triangles. In the analysis of the diffusion equation only triangulation with 4000 
and 16000 triangles are used because of the large measurement times. 
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2.1 Geodesic distance 



In the dynamical triangulation the volume V is identified with the number of 
triangles N, while the geodesic distance R is done with the number of dual 
links r of geodesic curve defined as the shortest path on dual link between two 
triangles. Introducing the lattice spacing parameter e, one can write 

V = Ne 2 , R = are 2 ", (21) 

where a is a dimensionless constant parameter. The parameter v is introduced 
again because X — R/V v is dimensionless. 

Let Sm{t) be the number of triangles on the boundaries which can be 
reached from a given triangle with the geodesic distance r. The average of 
Sn{t) under the dynamical triangulation, (Sn(t)), can be considered as a unit 
volume with the direction of r on fractal surface, and then satisfies 

oo 

X>iv(r)) = N. (22) 

The discrete version of U(X) is 

u(x) = N'-^SNir)) = x d »- l f(x) with x = ^ , (23) 

where /(0) is non-zero finite. u(x) and x are related with U(X) and X as 

u{x) = aU(X), ax = X. (24) 

By choosing a reasonable value of v, the same function u(x) is expected for 
any size of systems. Thus, one can determine the value of v from the point of 
view of finite size scaling. One can also determine the value of dh by 

db f M(x) = d h - 1 + const. x k + o(x k ) . (25) 
dlogx 

In refs. |0,tt3], <ilog (iSjv(r)) /dlogr for large N is studied by numerical 
simulation. In rcf. |]l2f the authors have succeeded in observing the scaling law 
for c = — 2 model by using quite a lot of triangles (five million triangles). We 
will show that the universal function u(x) helps us much more effectively to 
observe the fractal structure in numerical simulations. 
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Figure 1: u(x) as well as d log u(x)/d log x for pure gravity with system sizes IK, 2K, 4K, 

32K triangles as a function of the scaled variable x = (r + 5.5)/(Af 1 / 4 - 0.45). The 
theoretical distribution (the fully drawn line) is also displayed. 



Pure gravity 

Pure gravity is a good test case for numerical simulations because we know the 
exact formula for Gy (R) in ( 10 ) , the value of v = 1/4 and a = \Jq/{12 + 13\/3), 



and consequently the exact expression for u(x). Comparing the numerical data 
with the analytic curve, we find that a much better result is obtained by using 
a "phenomenological" variable x, 

(26) 



+ b ' 

instead of x = r/N 1 '. This shift is reasonable from the point of view of lattice 
artifacts. 

In fig. |j] we show the data as well as the theoretical curve of u(x) with an 
optimal choice of a and b. The agreement is almost perfect and the conclusion 
is that we already see continuum physics for systems as small as 1000 triangles 
in the case of pure gravity if we include simple finite size corrections like (|26|). 



Gravity coupled to matter 

We can now perform the analysis outlined above in the case of gravity coupled 
to matter: Ising spins coupled to gravity (c = 1/2), three-state Potts model 
coupled to gravity (c = 4/5) and one to five Gaussian fields coupled to gravity 
(c = 1, . . . , 5). The matter fields are placed in the centers of triangles. In these 
cases we do not know the theoretical values of v and dh ■ 

For all these theories we measure u(x) and try to determine the val- 
ues of v and dh according to (p3|), ( |25| ) and (p6|). We illustrate u{x) and 
dlogM(a;)/c?loga; for c = 1/2, 1, 2 and c = 5 in fig. H. Each of the graphs 
contains the scaled data for system sizes 1000, 2000, 4000, 8000, 16000 and 
32000 triangles. The results for other models (c = 4/5, 3 and 4) are similar. If 
v is chosen to be 1/4 for c < 1, 1/3 for c = 2 and 1/2 for c = 5, respectively, 



6 



we see as good scaling as in the case of pure gravity. For c < 1 the best values 
of the constants a and b in ( p6| ) are very close to the pure gravity values. The 
result of \ jv is shown on the left in fig. |||] From the right in fig. g we find that 
the "smooth" fractal condition (Q) is almost perfectly satisfied for all these 
theories. We have also checked that the exponent k defined in (|25| ) has the 
value, k ~ 4 for c < 1 while k ~ 2 for c = 4 and 5. 

We get the string susceptibility 7 S easily by the benefit of the minbu surgery 
algorithm. On the right in fig. we show the measured 7 S for various theories. 
Both figures in fig. ^| corroborates on the idea that 2D quantum gravity with 
large c is the branched polymer model which has \ jv = dh =2 and 7 S = 1/2. 

2.2 Time of Diffusion 

We now turn to the measurement of the spectral dimension d s introduced in 
( |l8| ) and ( |T^ ) . The diffusion equation at the discrete level describes a random 
walk. We denotes the discrete version of Kv(R) T) by fcAf(r; t). £ is the number 
of step of random walk, which is related with T as 

T = [3ts 2X , (27) 

where /3 is a dimensionless constant parameter. Since fcj\r(r;t) satisfies 

oo 

J2{SN(r))k N (r;t) = 1, (28) 

r=0 

the discrete version of ( |l6|) is 

- 1 r t 

kN(r;t) = —p(x;y) with x = — , y = —. (29) 

p(x;y) is related with P(X;Y) as 

p(x;y) = P(X;Y), [3y = Y . (30) 
The return probability of random walk is 

N \d B /2~l 

fcjv(0;i) ~ const. — ^d~/2 — ^ or * ^ , (31) 
while the average geodesic distance travel by random walk over time t is 

oo 

(r N (t)) d M ^(S' A r(r))rfc Ar (r;t) - const. N u ~ Xtx t a for i~0. (32) 



r=0 



d Closely related work has recently appeared in fllCfl. However, the situati on i n pure gravity 
is actually much better than it appears from the Taw" data presented in ||l(J. 
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Figure 2: u(x) as well as d log u(x)/d log x for c = 1/2, 1, 2 and 5. We have assumed 
v = 1/4 for c = 1/2 and c = 1, v = 1/3 for c = 2 and u = 1/2 for c = 5. The constants a, b 
determined this way are quite close to the pure gravity values for c = 1/2 and c = 1. In all 
cases the data include the following discretized volume sizes: N = IK, 2K, 4K, 32K. 
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Figure 3: The left figure shows l/u (as dh) determined by finite size scaling for c = 0, 1/2, 
4/5, 1, 2, 3, 4, and 5. The dots denote the best values of l/u d^). The right figure 
shows the measured string susceptibility ■js versus central charge c. We here assume 
Zy = const. V Jb ~ 3 . This is the reason why the measured 7 3 for c = 1 disagrees with the 
theoretical value. 




20 40 SO BO 100 20 40 60 SO 100 



Figure 4: The left figure shows — 2d log fcjv(0; t)/d log t versus t for c = (top curve), 
c = 1/2, c = 1, c = 3 and c = 5 (bottom curve). The system size is N = 16K. The right 
figure shows (r'jv(t)) 4 versus t for c = (bottom curve), c = 1, c = 3 and c = 5 (top curve) . 
The system size is N = 4K. Straight lines (as observed) indicate a = 1/4 according to ( p2j ) 

In fig. H we show — 2d\ogkN(0',t)/dlogt d s ) and (rAr(i)) 4 as a function of 
time for c = (pure gravity), 1/2, 1, 3, 5. From the left in fig. |I] we observe 
that d s is consistent with 2 for c < 1 and that it decreases for c > 1. From the 
right figure we observe that a « 1/4. We then obtain dh ~ 2d s if the "smooth" 
fractal condition ( po|) is satisfied. 



3 Discussion 

We generalize the definition of the two-point function in (||) as follows: Let 
^g{£i £o; Mi) be an observable which depends on two coordinates £ and £o 
as well as the metric, where Mj symbolizes some parameters. We define the 
average of ^f g (^, £o; Mi) with a fixed geodesic distance R as 
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x yy d 2 tVmd%yM^)5(d g (t, £ ) - R) Co; Mi) . 

Note that M f ) satisfies 

rdR{Sv(R))*v(R;Mi) = j (34) 

Jo 17 

where 



x //A\/9(£)<i 2 &\/9(&)*»({,&;M,). 

If (^(Mi))y 7^ 0, one can define a universal function, 

Q(X;Ft) = (*(M0)v * y(fl;M<) with X = v*' Yl = V^> (36) 

where a parameter /i, makes Yi = Mj/V*' dimensionless. Similarly to (|^) and 
@, we find 

\ dXU(X)Q(X;Y l ) = 1. (37) 
Jo 

The universal function of the moment like (|l9|) is in general expressed as 

o 

dXU(X)f(X;Y i )Q{X:Y l ) 1 (38) 







which is obtained by multiplying /(<f 9 (£, fo); M^) in front of ^f g in (p3|). It is 
straightforward to obtain the discrete version of ( pRf ) , ( p7| ) and (|||) . 

We have several examples as follows: i) = 1 is a trivial case because 



Q{X) — 1. ii) ^ g = K g (£, £o',T) is already discussed in section (1.2), i.e., 
Q{X;Y) = P(X;F). iii) * 9 = MOMHo) leads to Q(X) = V A ^+ A ^ V (R), 
where Aj is defined by {$\/g<j>i) v oc V 1 ~ Ai . If the surface is a "smooth" 
fractal, one can expect ^v(R) ~ const. iJ-(Ai+A a )/v for # _ . 

Acknowledgments 

It is a pleasure to thank Jan Ambj0rn and Jerzy Jurkiewicz for many interest- 
ing discussions. 

10 



References 



1. J. Ambj0rn and Y. Watabiki, Nucl. Phys. B 445, 129 (1995). 

2. J. Ambj0rn, J. Jurkiewicz and Y. Watabiki, Nucl. Phys. B 454, 313 
(1995). 

3. J. Ambj0rn, B. Durhuus and T. Jonsson, Phys. Lett. B 244, 403 (1990). 

4. Y. Watabiki, Prog. Theor. Phys. Suppl. No. 114, 1 (1993). 

5. H. Kawai, N. Kawamoto, T. Mogami and Y. Watabiki, Phys. Lett. B 
306, 19 (1993). 

6. J. Ambj0rn, P. Bialas, J. Jurkiewicz, Z. Burda and B. Petersson, Phys. 
Lett. B 325, 337 (1994). 

7. J. Ambj0rn and J. Jurkiewicz, Nucl. Phys. B 451, 643 (1995). 

8. S. Jain and S. Mathur, Phys. Lett. B 286, 239 (1992). 

9. J. Ambj0rn, S. Jain and G. Thorlcifsson, Phys. Lett. B 307, 34 (1993); 
J. Ambj0rn, G. Thorlcifsson, Phys. Lett. B 323, 7 (1994). 

10. S. Catterall, G. Thorlcifsson, M. Bowick and V. John, Phys. Lett. B 354, 
58 (1995). 

11. M.E. Agishtein and A. A. Migdal, Int. J. Mod. Phys. C 1, 165 (1990); 
Nucl. Phys. B 350, 690 (1991). 

12. N. Kawamoto, V.A. Kazakov, Y. Saeki and Y. Watabiki, Phys. Rev. Lett. 
68, 2113 (1992); Nucl. Phys. B (Proc. Suppl.) 26, 584 (1992). 



11 



